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O : ABSTRACT 

We have constructed a comprehensive grid of 690 metal line-blanketed, NLTE, 
plane-parallel, hydrostatic model atmospheres for the basic parameters appro- 
priate to O-type stars. The Ostar2002 grid considers 12 values of effective 
temperatures, 27 500 K < T eff < 55 000 K with 2 500K steps, 8 surface gravi- 
ties, 3.0 < log g < 4.75 with 0.25 dex steps, and 10 chemical compositions, from 
metal-rich relative to the Sun to metal-free. The lower limit of logg for a given 
effective temperature is set by an approximate location of the Eddington limit. 
The selected chemical compositions have been chosen to cover a number of typical 
environments of massive stars: the galactic center, the Magellanic Clouds, Blue 
Compact Dwarf galaxies like IZwl8, and galaxies at high redshifts. The paper 
contains a description of the Ostar2002 grid and some illustrative examples 
and comparisons. The complete Ostar2002 grid is available at our website at 
http : / /tlusty . gsf c . nasa . gov. 

Subject headings: stars: atmospheres, early-type — methods: numerical — radiative 
transfer 
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1. Introduction 

Theory of stellar atmospheres is presently enjoying a period of revival. This renewed 
interest has several major reasons. 

On the observational side, there has been an unprecedented expansion of quantity and 
quality of stellar spectrophotometric observations. Thanks to high sensitivity and high 
signal-to-noise ratio, one can now routinely record good-quality spectra of individual stars 
in other galaxies than the Milky Way, thus giving birth to a field of extragalactic stellar 
astronomy (Kudritzki 2000). The importance of deriving accurate stellar parameters as 
significant constraints for the global galactic and cosmological parameters, coupled with the 
quality of observations, compels us to use as accurate and reliable model atmospheres as 
possible. 

Among all stellar types, the hot, massive, and bright stars - in particular O-type stars, 
stand as one of the most important stellar groups in the global astronomical context. They 
contribute most of the radiation to the UV and optical integrated spectra of galaxies which 
still undergo a star-forming process. The O stars provide the majority of ionizing photons 
for H II regions, and the interstellar and intergalactic medium. In the cosmological context, 
there is still a debate whether hot massive stars or Active Galactic Nuclei provided the 
necessary photons for the re-ionizing of the Universe around z ~ 3. In view of forthcoming 
space missions, like the New Generation Space Telescope, there is a growing interest in 
studying the early (or even the first) generation of stars in the Universe (see, e.g., Weiss 
et al. 2000). 

On the modeling side, the capabilities of the modeling methodology has advanced enor- 
mously. One obvious reason is a many-fold increase of computer speed and memory. Another, 
and also very important, reason is the development of efficient and robust numerical meth- 
ods. The central role in this respect is played by the Accelerated Lambda Iteration (ALI) 
method (for a recent review, see Hubeny 2003). 

With the development of fast computers and efficient numerical methods, the basic 
limiting factor is now becoming the accuracy (and availability) of atomic data. Fortunately, 
great advances are being made in this direction too. Recently, two major collaborative 
projects, the Opacity Project (OP) (Seaton 1987; The Opacity Project Team 1995; 1997), 
its continuation the IRON Project (Hummer et al. 1993; Pradhan et al. 1996; Nahar 2003), 
and the OPAL project (Iglesias & Rogers 1991, 1996), have produced atomic data on a large 
scale. 

The time is thus ripe to compute a new grid of model stellar atmospheres. We need to 
analyze spectra of individual stars as well as composite spectra of galaxies, with much larger 
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degree of complexity and accuracy than is provided by the Kurucz (1993) spectral energy 
distributions (SED's) currently in use. 

Because of the importance of O-type stars, such a modeling effort is hardly new. How- 
ever, there are three issues that make a construction of O star model atmospheres a highly 
non-trivial task. They are (i) the presence of a strong mass outflow, the stellar wind; (ii) 
significant departures from Local Thermodynamic Equilibrium (LTE); and (Hi) an effect of 
numerous metal lines, traditionally called metal line blanketing. At present, there are no 
models that would address all these three issues in full (that is to construct fully blanketed, 
non-LTE models for an entire atmosphere, going from essentially static photosphere out to 
the wind - the so-called unified models), although there are several research groups and 
several computer programs that are in principle capable of doing so - programs CMFGEN 
(Hillier & Miller 1998; Hillier 2003); Hoeflich's suite of programs (Hoeflich 1995, 2003); 
PHOENIX (Hauschildt et al. 1997); Munich programs (Pauldrach et al. 2001); and the 
Kiel- Potsdam code (Hamann 1985; Koesterke et al. 2002). 

As we will discuss in more detail in the next Section, unified models either require huge 
computer resources, or one has to approximate the other two ingredients - NLTE and line 
blanketing. Also, a description of a stellar wind is still somewhat uncertain. It is therefore 
desirable to be able to construct complete models for hydrostatic layers, traditionally called 
the stellar photosphere. Even this simplified problem is a tremendous task. Early NLTE 
model atmospheres (Auer & Mihalas 1969b; Mihalas 1972) considered only a few, about 10, 
energy levels of H and He. Likewise, only most important opacity sources were considered - 
several bound-free transition of H and He, and a few most important hydrogen lines. On the 
other hand, both observations and LTE line-blanketed models (Kurucz 1979) showed that the 
UV spectrum contain a huge number of metal, mostly iron, lines, which may then significantly 
influence not only the emergent spectral energy distribution, but the atmospheric structure 
in general. 

Although it was clear that early NLTE models provided a significantly improved pre- 
dictions of line profiles of most important spectral lines (see, e.g., a number of papers in 
Garmany 1990; Crivellari, Hubeny, & Hummer 1991; Heber & Jeffery 1992) it was not clear 
what is most significant from the point of view of an overall model atmosphere accuracy - 
to compute NLTE models without influence of metal line blanketing (or with a very simpli- 
fied blanketing), or to treat metal line blanketing as completely as possible, while making a 
sacrifice of neglecting NLTE effects. 

Anderson (1985) has introduced an ingenious multi-frequency/multi-grey method, and 
was able to construct models with several tens to hundreds of energy levels and hundreds 
of lines treated explicitly. When the ALI method was first implemented (Werner 1986), 
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the number of lines that can be treated explicitly increased to thousands. Anderson (1989) 
used a concept of superlevels and superlines to describe effects of numerous lines of iron- 
peak elements. Later, even more complex model atmospheres with improved treatment of 
metal line blanketing were constructed (Dreizler & Werner 1993; Hubeny & Lanz 1995). The 
experience gained during the last decade taught us that both NLTE and line blanketing are 
very important; moreover, one needs to treat both to a very high degree of completeness in 
order to avoid spurious results. 

With the recent improvements of numerical schemes and with extensive developments 
of computational resources that have occurred during the last decade we are now able to 
compute essentially "exact" , fully blanketed NLTE hydrostatic model atmospheres. To fulfill 
this goal is the topic of the present paper. 

The paper is organized as follows. In § 2 we summarize the basic physical assumptions 
and the numerical methods used to construct our model grid, and in § 3 we specifically 
describe numerical approaches to treat the metal line blanketing. Section 4 is devoted to an 
overview of atomic data, and, in some cases, their approximations. Section 5 describes the 
Ostar.2002 model grid in more detail. Some representative results are presented in § 6, while 
global parameters such as the number of ionizing photons and a bolometric correction are 
discussed in § 7. Section 8 discusses a sensitivity of computed models to adopted parameters. 
Our results are summarized in §9. 

2. Assumptions and Methods 

Our basic physical assumptions are those of plane-parallel stratification, hydrostatic 
equilibrium, and radiative equilibrium. As has been pointed out in the literature repeatedly 
over the past three decades (Auer & Mihalas 1969b; Mihalas 1978, a number of papers in 
Garmany 1990; Crivellari, Hubeny, & Hummer 1991; Heber & Jeffery 1992; and Hubeny, 
Mihalas, & Werner 2003), the approximation of Local Thermodynamic Equilibrium (LTE), 
which is another typical approximations of classical model atmospheres, breaks down for 
early-type stars, in particular for O stars. We thus drop the approximation of LTE and 
compute non-LTE (NLTE) models. 

2.1. Are hydrostatic O-star model atmosphere relevant? 

At the very outset, we have to clarify one crucial question. It is well known that the most 
common features of the UV spectra of O stars are strong P Cygni profiles of numerous metal 
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lines (C IV, N V, O V, Si IV, S IV, S V, P V, etc.), which provide an unmistakable evidence 
of massive, supersonic outflows - the stellar winds. Since our models are hydrostatic, we 
have to ask: are the hydrostatic model atmospheres for O stars relevant at all? 

The answer depends on what exactly we are expecting from a model. If we intend to 
analyze the P Cygni profiles, and to derive basic wind parameters (mass loss rate, terminal 
velocity), then, obviously, the hydrostatic models at best provide lower boundary conditions. 
However, if we intend to derive the basic stellar parameters (effective temperature, surface 
gravity - therefore mass - , and chemical composition, from the bulk of the UV and optical 
spectrum), then the answer is that the hydrostatic models are not only useful, but in many 
cases actually preferable over the wind models or the unified (photosphere + wind) models 
with a simplified treatment of a photosphere. There are several reasons for this assessment. 

First, most spectral lines (excluding strong, mostly resonance, lines) are formed in the 
photosphere where velocities are small and geometrical extension is negligible. Indeed, the 
spectrum of slowly-rotating O stars show mostly narrow and symmetric lines (see, e.g., the 
high-resolution UV spectral atlas of the 09V star, 10 Lac, Brandt et al. 1998; and the 
quality of the fit of the HST spectra by our NLTE line-blanketed models - Hubeny, Heap, & 
Lanz 1998). This statement is valid in all but the most extreme cases of stellar winds, like 
those observed in Wolf-Rayet stars and extreme Of supergiants. Moreover, winds are weaker 
in low-met allicity environments making photospheric models even more relevant in this case. 
The theoretical spectra predicted by the wind code CMFGEN (Hillier & Miller 1998) and 
by Tlusty agree very well for a test case, T c g=35 000 K, log g = 4.0, demonstrating that 
most lines are indeed formed in the quasi-static photosphere (Hillier & Lanz 2001; Bouret 
et al. 2003). 

Second, there are many remaining uncertainties as to the exact properties of stellar 
winds of hot stars, even though the paradigm of radiatively-driven winds is well established. 
For example, we do not have in many cases a consistent solution: the calculated line force 
is often too small to drive a wind through the critical point (though some progress is being 
made), so a /3-type velocity law must be adopted to describe empirically the velocity and 
density structure of the wind. Furthermore, radiative equilibrium is not really satisfied 
in these winds, where shocks dump mechanical energy and heat the wind, thus changing 
its ionization structure (e.g., super-ionization). Finally, the assumption of one-dimensional 
geometry is challenged by the likely presence of dense clumps in the wind and by rotation; 
the role of magnetic fields is not yet understood, and their importance in structuring the 
wind will very likely be more important than in the photosphere. 

Third, there is also a practical aspect. Spherical model atmospheres with velocity 
fields, although technically still a 1-D problem, have effectively a larger dimensionality than 
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a plane-parallel static problem because of a sharp angular dependence of the radiation field. 
(Recall that in spherical wind models, the number of angles is about the same as the number 
of discretized radial points, while in the plane-parallel case considering a few - typically 
3 - angles is usually sufficient). Consequently, computing spherical expanding models is 
intrinsically more complex and more computationally intensive than plane-parallel models. 
One is thus forced either to compute fewer models, or to make approximations in the local 
physics (in particular opacities), to reduce computer requirements to reasonable limits. 

Finally, it is important to realize that while line blanketing operates everywhere in the 
atmosphere, departures from hydrostatic equilibrium occur only in the outer layers. His- 
tory has taught us that it is very dangerous to neglect or even to approximate metal line 
blanketing. Neglecting line blanketing leads not only to a lower local temperature in the 
continuum-forming layers and thus to systematic errors in determining the effective temper- 
ature (e.g. Hubeny, Heap & Lanz 1998), but also to errors in predicted ionization balance of 
important elements (CNO, etc.). These species provide in turn important diagnostic tools 
for the wind. Since the ionization edges of CNO and other light metals are located in the 
far and extreme UV, and since this region is crowded with numerous lines of highly ionized 
iron and nickel (Fe IV- Fe VI; Ni IV-Ni VI), any deficiencies in describing photospheric line 
blanketing are directly reflected in the uncertainties in the predicted wind structure. It is 
thus no exaggeration to say that without a proper and detailed treatment of photospheric 
metal line blanketing even sophisticated wind models may be seriously in error. 

Therefore, we are on much safer ground to determine the fundamental parameters of 
massive stars as long as we use selected photospheric lines and avoid all uncertainties re- 
maining with wind models. Moreover, a grid of NLTE line-blanketed model atmospheres 
incorporating state-of-the-art opacities will provide a sound basis for comparison to future, 
more realistic model atmospheres of O stars. 

2.2. Numerical Methods 

Over the years, we have developed a computer program called Tlusty for computing 
model stellar atmospheres without assuming LTE (Hubeny 1988; Hubeny & Lanz 1992, 
1995, 2003; Lanz & Hubeny 2001; Lanz et al. 2003). The program solves the basic equations 
- radiative transfer, hydrostatic equilibrium, radiative equilibrium, statistical equilibrium, 
charge and particle conservation. 

These equations are discretized in frequency and depth, which yields a set of highly- 
coupled, non-linear algebraic equations. We have used in this work the so-called hybrid 
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CL/ALI method that was introduced by Hubeny & Lanz (1995). This method builds on 
the advantages of both ALI and Complete Linearization. Through ALI, the mean intensities 
of the radiation field are eliminated from the set of linearized equations. However, the hy- 
brid CL/ALI method still allows us to linearize the radiation intensity at a few selected key 
frequencies, notably improving the convergence rate compared to a "pure" ALI approach. 
The hybrid CL/ALI method combines therefore the low cost of ALI while retaining a con- 
vergence rate that is similar to CL. Moreover, levels and superlevels may be grouped during 
the linearization step, reducing even further the number of linearized equations (Hubeny & 
Lanz 2003). The full statistical equilibrium is solved in a separate step, element by element, 
with a typical maximum number of levels around 200 per element. 

Under the above assumptions, the critical ingredient of the models is a degree of com- 
pleteness of the opacity sources, and of the atomic/ionic energy levels treated explicitly in 
the statistical equilibrium equations. We shall discuss these issues in more detail in the next 
two Sections. 



3. Treatment of Metal Line Blanketing 

As we have pointed out above, NLTE line blanketing is an essential ingredient of real- 
istic model atmospheres, determining the vertical temperature structure and the ionization 
structure. NLTE ionization rates are typically larger than their LTE counterparts, because 
the radiative rates are determined by the radiation field which is emitted in deeper, hotter 
layers. Therefore, accurate mean intensities of the radiation field are required to evaluate 
correctly the NLTE ionization structure, which in turn require an accounting of all sources 
of opacity. When neglecting in particular the opacity of Fe lines, NLTE models predict 
too strong EUV/UV radiation fields, and overestimate the ionization and the importance 
of departures from LTE. In order to construct realistic NLTE model atmospheres, we have 
therefore included about 100,000 individual atomic levels of over 40 ions of H, He, C, N, O, 
Ne, Si, P, S, Fe, and Ni in the calculations. The levels are grouped into about 900 NLTE 
superlevels, with populations of individual energy levels following Boltzmann statistics in- 
side a superlevel. All individual levels in a superlevel thus share the same NLTE departure 
coefficient. Hubeny & Lanz (1995) discussed the concept of NLTE superlevels in detail. A 
total of 8000 lines of the light elements and about 2 millions lines of Fe m-vi and Ni ill-vi 
are accounted for in the calculations. Based on a line strength criterion, Fe and Ni lines are 
dynamically selected from a total of 12 million lines listed by Kurucz (1994). The numbers 
of NLTE levels and lines included in the model atmospheres are listed in Table 1. Details of 
the model atoms are given in Sect. 4. 
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ALI and level grouping are two major steps toward the practical realization of NLTE 
line-blanketed model atmospheres. "Superlines" , i.e. transitions between superlevels, are 
another important ingredient. A superline may involve hundreds or thousands of individual 
lines over a relatively broad frequency range. The resulting absorption cross-section is then 
very complex and a detailed representation may require a very large number of frequency 
points. In an early implementation of NLTE line blanketing, Hubeny & Lanz (1995) adopted 
the concept of Opacity Distribution Functions (ODFs) to describe the absorption profile of 
superlines. The total opacity of all lines in a given transition (computed using Voigt profiles 
for each line) is sorted; the sorted opacities could then be represented with a limited number 
of frequencies, typically 15 to 30 per transition, or 30 000 to 50 000 frequencies for the 
whole spectrum. The ODF approach assumes that the background opacity changes little 
over the frequency range covered by the superline. Details of blends with other lines or 
continua are therefore lost; however, this issue should not be serious in a statistical sense 
since this method is designed for millions of lines of the iron-peak elements. One may still 
unintentionally create or neglect blends between strong lines and ODFs, thus affecting the 
resulting model atmospheres. 

We have therefore implemented a second method, Opacity Sampling (OS), that have 
the advantage of treating exactly the blends. OS is a simple Monte Carlo-like sampling 
of the superline cross-sections. Contrary to the ODF approach, OS may miss the cores of 
important strong lines if too few frequency points are used. We use a variant of OS, sampling 
the whole spectrum at prescribed intervals in frequency. The typical sampling step is 0.75 
fiducial Doppler widths of Fe. Smaller steps are occasionally used to include light elements' 
line cores in the frequency set as well as points around photoionization edges. Our model 
atmospheres are sampled with 180 000 to 200 000 frequency points for an O star model. 
While costly in terms of computing requirements, we are thus sampling the spectrum with 
sufficient resolution to ensure that every line is accounted for. This is the current standard 
mode in Tlusty, which was used to compute this grid of model atmospheres. 

Lanz & Hubeny (2003) presented some illustrations of the ODF and OS approaches. 
In Sect. 8, we show the changes in a typical model atmosphere that result from different 
representations of line opacities. 

4. Atomic Data 

The Opacity Project (OP 1995, 1997) provides very extensive datasets for radiative 
transitions (^/-values, photoionization cross-sections) which have been obtained from ab- 
initio calculations for all ions of the most abundant light species (Z < 14, S, Ca, and Fe). 
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The current status of the project and its continuation as the Iron Project is discussed by 
Nahar (2003). OP claims that their data are accurate at the 10% level or better. We have 
adopted systematically these data for light elements and they form the bulk of radiative 
data used in Tlusty. The OP data have been extracted from Topbase 2 , the database of 
the Opacity Project (Cunto et al. 1993). OP also provides theoretical level energies. Despite 
the sophistication of OP calculations, these energies still differ by a few percent from the 
energies measured in the laboratory. We have therefore updated the level energies with 
experimental energies when they are available in the Atomic and Spectroscopic Database 3 
at NIST (Martin et al. 1999). These large datasets are manipulated with Modion, an IDL- 
based graphic tool that works directly with Topbase data files (Varosi et al. 1995), to prepare 
the model atoms required by Tlusty. For most lines, we have adopted depth-independent 
Doppler line profiles, with a standard temperature, T=0.75T o g. Depth-dependent, Stark 
or Voigt profiles are used for the strongest lines (about 70 lines of light elements and the 
lower line series of H I and He II) . Radiative damping constants are extracted from the NIST 
database, and Stark widths are adopted from the work of Dimitrijevic and collaborators. 
OP photoionization cross-sections contain many autoionization resonances. As Tlusty 
handles large frequency sets, we could therefore adopt the theoretical OP cross-sections 
in every detail. However, sharp resonances are not fully resolved by the OP calculations. 
Moreover, resonances are shifted off their exact wavelengths due to the limited accuracy of 
theoretical energies. Therefore, we have followed Bautista et al. (1998) and implemented the 
concept of Resonance-Averaged Photoionization (RAP) cross-sections. The cross-sections 
are smoothed by a Gaussian with a width, 5E/E = 0.03, that corresponds to the typical 
uncertainty of the resonance energies. OP photoionization cross-sections are total cross- 
sections for each initial state summed over all final states. We have kept this approximation: 
we assume that all bound-free transitions considered by Tlusty end up in to the ground 
state of the next ion. If the transition to the ground state is forbidden by selection rules, 
the cross-sections is shifted in energy corresponding to the excitation of the lowest excited 
state for which the transition is allowed. While the opacity and the ionization rate are 
thus computed correctly, the corresponding recombination rate (based on detailed balance 
arguments) is only approximate because it assumes that the end state of ionization is in 
Boltzmann equilibrium with the ground state. Nevertheless, since the population balance of 
the upper ion is determined by a large number of processes (mostly bound-bound transitions), 
the influence of an approximate recombination rate is typically small, and the resulting error 
in the statistical equilibrium virtually negligible. Additional details on the model atoms for 



2 http : //heasarc . gsf c . nasa . gov/topbase/ 

3 http : //physics . nist . gov/cgi-bin/AtData/maimisd 
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each ion is given in the relevant subsections. 

For iron-peak elements, we have extracted level energies and line oscillator strengths 
from Kurucz (1994) extensive semi-empirical calculations. Kurucz (1994) lists thousands of 
energy levels for each ion. A statistical approach to deal with so many levels is therefore 
necessary; solving the full statistical equilibrium equations would indeed require inverting 
very large matrices which is not only costly but numerically inaccurate. We therefore group 
individual levels into superlevels, assuming that the population of each level in a given 
superlevel follow the Boltzmann distribution. While the populations of superlevels could 
depart from their LTE values, all sublevels in a given superlevel share the same 6-factor. We 
have applied two criteria to group the levels: all sublevels must have very close excitation 
energies, and they share the same parity. The first criterion ensures that collisional rates 
between sublevels are large enough to achieve the Boltzmann statistics; the second criterion 
further supports the validity of this assumption by treating properly transitions out of and in 
to metastable states. The second requirement provides the additional advantage of avoiding 
radiative transitions within a single superlevel. The absorption cross-section of a superfine 
between two superlevels can be quite complex, involving hundreds or thousands of individual 
lines. Exact absorption cross-sections are computed at three depths (top and bottom of the 
model atmosphere, and at a depth where Tr oss ~ 1 and T ^T cS ), and are logarithmically 
interpolated between these points. A detailed representation requires a very large number 
of frequency points. Opacity Sampling is a method designed for handling such cases, and its 
implementation in Tlusty was outlined in Sect. 3. 

Non-LTE calculations also require collisional excitation and ionization rates. However, 
collisional rates are known only for a limited set of transitions, mostly inter-system transi- 
tions that are important for diagnostics at low densities. The Iron Project (Hummer et al. 
1993) is starting to remedy this situation, but their data are not yet implemented in Tlusty. 
We rely on general formulae at this stage. We use the Seaton (1962) formula for collisional 
ionization; for optically- allowed line transitions, we use the van Regemorter (1962) formula 
which is accurate to about 30 - 50 % (Pradhan, priv. comm.); for forbidden transitions, no 
approximate formula or a simple prescription is available, so we somewhat arbitrarily use 
Eissner-Seaton's formula with 7(T) = 0.05. Different approaches of collisional strengths in 
superlines may be implemented. Lanz & Hubeny (2003) compared two different estimates 
and they argued that the most natural way consists in summing up all contributions for indi- 
vidual transitions, using the appropriate formulae for allowed and for forbidden transitions. 
We have adopted this approach here. 

The following subsections provide the relevant details of the model atoms. Table 1 
summarizes the atomic data included in the model atmospheres, as well as the references 



- 11 - 



to the original calculations collected in Topbase. Datafiles and Grotrian diagrams may be 
retrieved from TLUSTY's Web site, http://tlusty.gsfc.nasa.gov. 

4.1. Hydrogen 

The model atom of hydrogen includes the lowest 8 levels and one superlevel grouping 
the higher excitation levels, from n = 9 to n = 80, following the occupation probability 
formalism developed by Hummer & Mihalas (1988) and extended to NLTE situations by 
Hubeny et al. (1994). The Lyman and Balmer continuum cross-sections are generalized to 
account for transitions from these two levels to dissolved higher states. This contribution is 
neglected for all other continua. 

All lines between the first 8 bound levels are included, assuming an approximate Stark 
profile for the Lyman and the Balmer series (Hubeny et al. 1994) and a Doppler profile for the 
other lines. The higher members of the Lyman and the Balmer line series are incorporated 
via the OS approach. This is the only difference from the treatment of Hubeny et al. (1994) 
who adopted ODFs. OS is used here in order to account exactly for line blending with the 
He II line series (and with other lines). Higher members of other line series are neglected. 

Hydrogen is the main opacity source in the photospheres of O stars. We can benefit 
most of our hybrid CL/ALI method by selecting a few key frequencies where the hydrogen 
opacity is largest, at which the mean intensity of the radiation field is linearized. We have 
selected the first six frequencies immediately higher than the Lyman edge, the three core 
frequencies in Ly a, and the core frequency of Ly j3, Ly 7, H a, and H (3. A few more 
frequencies have been selected in the cores of strong resonance lines of other species, and 
will be mentioned in the relevant subsections. Although there might be a different optimal 
choice of these frequencies for each individual model atmosphere, we have found that this 
choice was generally adequate to achieve most of the gains of the CL/ALI method in the 
parameter space covered by the Ostar2002 grid. 

4.2. Helium 

Neutral helium is modelled with a 24-level atom, which incorporates all 19 individual lev- 
els up to n = 4, and 5 superlevels grouping levels with identical principal quantum numbers 
(5 < n < 8); for n = 5, singlet and triplet levels are grouped into two separate superlevels. 
We have performed some tests with more detailed model atoms, splitting the superlevels up 
to n = 6, or including higher excitation levels, but we did not find significant differences. 
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All lines from n < 4 levels are included, but infrared lines between high-excitation levels 
(n > 5) are assumed to be in detailed radiative balance. Voigt line profiles are used for 
lower-excitation lines (n < 2), and Doppler profiles are assumed for higher-excitation lines. 
He I is treated in a special way in TLUSTY; like hydrogenic ions, and contrary to all other 
ions, we have special subroutines in the code to deal with atomic data. We use Hummer's 
(priv. comm.) routine to evaluate the photoionization cross-sections; the routine is based 
on Seaton & Fernley's cubic fits of OP data, and properly manages the different ways of 
averaging levels. Bound-free and bound-bound collisional transitions rates are also evalu- 
ated with a special routine based on approximate expressions from Mihalas et al. (1975). 
For bound-bound transitions between low-excitation levels, these approximate rates are su- 
perseded by detailed calculations of collisional cross-sections (Berrington & Kingston 1987), 
used by Storey and Hummer (priv. comm.) to evaluate collisional rate coefficients. 

The ionized helium model atom includes the first 20 bound levels. The He II model 
is built following closely H I. We will only indicate here the differences. The occupation 
probability formalism (Hubeny et al. 1994) is taking into account for the 20 bound levels, 
but higher-excitation levels are not included explicitly. We also use continuum cross-sections 
for bound-free transitions from the ground and first excited levels that account for transitions 
to highly-excited, dissolved levels. Stark profiles are adopted for the first four line series, 
while infrared lines (n > 10) are set to detailed radiative balance. Finally, the first three 
frequencies at the head of the He II Lyman continuum are linearized in the ALI/CL scheme. 



4.3. Carbon 

The general handling of atomic data has been discussed at the top of Sect. 4. We restrict 
the following subsections to a few particulars relative to each ion. Carbon is represented by 
a 22-level C II model atom, a 23-level C III model atom, and a 25-level C IV model atom, 
grouping 44, 55, and 55 individual levels, respectively. The C II model atom includes the 
first 17 levels, 14 from the doublet system and the first three levels of the quartet system. 
Five superlevels group higher-excitation levels, four in the doublet system (6 < n < 9) and 
one superlevel grouping three levels in the quartet system. The C III model atom includes 
the first 16 levels (E < 312 000 cm -1 ), and group all higher-excitation levels up to n = 6 
into 7 superlevels. Levels with n > 7 are not included explicitly (see further discussion in 
Sect. 8 about this restriction and its consequence on C 3+ recombination). The C IV model 
atom includes all individual levels up to n = 6 and four superlevels (7 < n < 10). The 
data extracted from Topbase have been extended to include all levels up to 10m 2 M°. We 
have included the two fine-structure levels of 2p 2 P° for an exact treatment of the C IVA1549 
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resonance doublet. Depth-dependent Voigt profiles are adopted for the resonance lines of 
the three ions, and for the strong C III A1176 line from the 3 P° metastable level. 

4.4. Nitrogen 

The model atmospheres incorporate four ions of nitrogen, N II to N V, and the ground 
state of N VI. The first 16 levels in the singlet and triplet systems of N II (n — 2; 3s, and 
3p levels), and the first two quintet levels, are included in the models as explicit individual 
NLTE levels. Higher-excitation levels (3d; 4 < n < 5) are grouped into 8 superlevels, four in 
each system. Higher levels (n > 6) are not incorporated explicitly in the model atom. The 
N III model atom includes all the levels below the ionization limit as extracted from Topbase. 
The first 24 levels (E < 330 000 cm" 1 ) are treated as individual levels, and the other 43 levels 
are grouped into 7 superlevels. The ground state is split into two fine-structure levels for 
an exact treatment of the resonance doublets (e.g., N IIIA990). The N IV model atom was 
constructed similarly to the C III model atom. The first 15 levels are included explicitly, and 
higher levels up to n = 6 are grouped into 8 superlevels. Levels with n > 7 are not included 
explicitly. Finally, the N V model atom is similar to C IV, although levels are grouped into 
6 superlevels starting with n = 5 levels. The two fine-structure levels of 2p 2 P° are included 
in the model atom. Depth-dependent Voigt profiles have been set up for 18 resonance lines 
of the four ions. 



4.5. Oxygen 

The model atmospheres incorporate five ions of oxygen, O II to O VI, and the ground 
state of O VII. O II and O III have a relatively rich structure of levels. High excitation 
levels have thus not been incorporated in the model atoms. The O II model atom includes 
21 individual levels and 8 superlevels, grouping levels up to n = 5, while the O III model 
atom includes the lowest 20 levels and groups higher-excitation levels up to n = 6 into 9 
superlevels. All O IV below the ionization limit are incorporated, the first 30 as individual 
levels, the others grouped into 8 superlevels. Fine-structure of the O IV ground state is 
included. The first 34 O V levels are included explicitly, and higher levels up to n = 7 are 
grouped into 6 superlevels. Finally, the O VI model atom is similar to C IV, although levels 
are grouped into 5 superlevels starting with n = 6 levels. The two fine-structure levels of 
2p 2 P° are included in the model atom. Depth-dependent Voigt profiles have been set up for 
O IVA789, 609, 554, and O VA630. 
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4.6. Neon 

Simple model atoms have been built for the neon ions, Ne II to Ne IV. We are not aiming 
at a detailed NLTE calculation of Ne populations, but our purpose was to account for the 
bound-free opacity in the extreme-ultraviolet due to the ionization of this abundant species. 
The Ne II model atom includes the first 11 levels and groups levels up to E < 300 000 cur 1 
into four superlevels. Ne III levels up to E < 400 000 cm -1 are included in the model atom, 
as 12 individual levels and 2 superlevels. Finally, the Ne IV model atom includes the first 
10 levels, and 2 superlevels grouping levels up to E < 550 000 cm -1 . All lines are assumed 
to have depth-independent Doppler profiles. 

4.7. Silicon 

Detailed model atoms have been constructed for Si III and Si IV, because the silicon 
ionization balance is a good temperature indicator for late O and B stars. The Si III includes 
all 105 energy levels listed by Topbase below the ionization limit. The first 24 levels are 
treated as individual explicit levels, while the higher levels are grouped into 6 superlevels 
(3 superlevels in the singlet and in the triplet systems). The Si IV model atom includes 
all individual levels up to n = 6 and four superlevels (7 < n < 10). The data extracted 
from Topbase have been extended to include all levels up to 10m 2 M°. We have included 
the two fine-structure levels of 2p 2 P° for an exact treatment of the Si IVA1397 resonance 
doublet. Depth-dependent Voigt profiles are adopted for Si III A1206, 567, and Si IVA1394, 
1403 resonance lines. 



4.8. Phosphorus 

Phosphorus has not been considered by OP, because the relatively low phosphorus 
abundance in the Sun indicates that phosphorus is most likely only a minor contributor to 
the total opacity in astrophysical plasmas. However, we have included phosphorus as an 
explicit NLTE species in our model atmospheres due to the potential diagnostic value of the 
P IVA951 and P V All 18-28 resonance lines. We have built simple model atoms based on 
data extracted from the NIST Atomic and Spectroscopic Database. The P IV model atom 
includes the first 14 levels, while the P V model atom includes individual levels up to n = 5 
and groups 6 < n < 9 levels into 4 superlevels. We have approximated the photoionization 
cross- sect ions by interpolating in cross-sections from the same configurations in isoelectronic 
sequences (e.g., using Si IV and S VI cross-sections to estimate P V cross-sections). We 
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believe that this approach is likely more realistic than a simple hydrogenic approximation. 

4.9. Sulfur 

The model atmospheres incorporate four ions of sulfur, S III to S VI, and the ground 
state of S VII. The model atoms are kept relatively simple, aiming at including the sulfur 
bound-free opacity in the model atmospheres and at predicting a few key diagnostics lines, 
mostly S IVA1073 and S VA1502. The S III model atom includes all levels with E < 
200 000 cm -1 , the first 16 levels individually and the other 11 levels grouped into 4 superlevels. 
The S IV model atom includes the first 14 levels up to E < 225 000 cm -1 , while the S V 
model atom includes the first 12 levels. The S VI model atom includes all individual levels 
up to n = 5 and three superlevels (6 < n < 8). Fine structure of S IV and S VI 3p 2 P° levels 
is taken into account. The S VIA933-45 resonance lines are the only sulfur lines for which 
we have adopted depth-dependent Voigt profiles. 

4.10. Iron 

The model atoms of iron-peak elements are based on Kurucz (1994) extensive semi- 
empirical calculations. We have extracted a list of even and of odd-parity levels for each ion 
included in the model atmospheres. We have considered all levels, those with experimental 
energies as well as those with predicted energies. We plotted histograms of numbers of levels 
per energy interval to set up the energy bands defining the superlevels. We chose typically 
15 to 30 bands per parity, so the model atoms have 30 to 50 superlevels. For each ion, we 
read the full linelist from Kurucz (1994) and select dynamically the strongest lines. The line 
selection criteria include line g /-values, excitation energies, and ionization fractions. In this 
way, we skip most or all Fe VI lines in the coolest models and most or all Fe III lines in the 
hottest models. We have extracted photoionization cross-sections from the latest calculations 
of Fe III to Fe V by the Ohio State group (see Table 1). They assumed LS'-coupling, and 
we could typically assign theoretical cross-sections to observed levels for the lowest 20 to 30 
levels. We have assumed an hydrogenic approximation for higher-excitation levels and for 
all Fe VI levels. The data were then summed up to setup cross-sections for the superlevels. 
Finally, we have applied to these cross-sections the RAP smoothing technique. 
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4.11. Nickel 

We have followed the same approach for nickel, using Kurucz (1994) data. OP did not 
include nickel in their calculations. We have therefore adopted an hydrogenic approximation 
for the photoionization cross-sections. 

5. Description of the Grid 

The Ostar2002 grid covers the parameter space of O stars in a dense and compre- 
hensive way. We have selected 12 effective temperatures, 27500 K<T c fj< 55 000 K, with 
2 500K steps, 8 surface gravities, 3.0 < log<? < 4.75, with 0.25 dex steps, and 10 chemical 
compositions, from metal-rich relative to the Sun to metal-free. Table 2 lists the selected 
metallicities. Solar abundances refer to Grevesse & Sauval (1998). In every case, we have 
assumed a solar helium abundance, He/H=0.1 by number. All other chemical abundances 
are scaled from the solar values. Tlusty can readily accommodate non solar-scaled abun- 
dances resulting, e.g., from mixing of CNO-cycle processed material; however, this would 
increase the total number of models to compute to an unmanageable level by increasing the 
dimensionality of the parameter space. Thus, we view this grid as a starting point from 
which detailed abundance analyses can be performed. The selected chemical compositions 
have been chosen to cover a number of typical environments of massive stars: the galactic 
center ("C models"), the Magellanic Clouds ("L models" for LMC; "S models for SMC"), 
Blue Compact Dwarf galaxies (e.g. "W models" for IZwl8), and galaxies at high redshifts 
("X and Y models"). 

For each composition, we have computed 69 model atmospheres, as illustrated in Fig. 1. 
Evolutionary tracks from the Geneva group (Schaller et al. 1992; Charbonnel et al. 1993) 
have been overplotted in Fig. 1. Stars with masses larger than about 15 M are in the range 
covered by our model atmospheres. The highest gravity and the lowest temperature have 
been included in the grid to avoid extrapolations when using this grid as a starting point for 
NLTE spectral analyses. 

The lower limit in log g is set by the Eddington limit. Using NLTE H-He model at- 
mospheres, we have estimated the lowest surface gravity as a function of T c s which would 
correspond to the Eddington limit neglecting the effect of rotation. This limit is reached for 
values typically 0.1 to 0.2 dex smaller than the grid limit. We present in Fig. 2 the variation 
of the radiative acceleration with depth for a solar composition model, T e ^= 40 000 K, and 
logg = 4.0. The radiative acceleration goes through a maximum in the continuum forming 
region (tr oss 1) and shows a strong increase in very superficial layers (tr oss < 10~ 6 ). Near 
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the surface, this model atmosphere becomes unstable (g Ta d > g)- Our modeling assumptions 
obviously break down in these layers. We have therefore numerically limited the radiative 
acceleration in superficial layers to ensure convergence. This case is quite typical; a similar 
behavior is observed in all our model atmospheres, although higher gravity models remain 
numerically stable at every depth (about two thirds of all solar composition models). The 
model spectra are not affected by this approximation, because these superficial layers only 
influence the strong resonance lines that form in a stellar wind. The stability of our model 
atmospheres is defined by the maximum reached by the radiative acceleration in the photo- 
sphere (see Fig. 2). Fig. 3 displays contours of this maximum relative to the gravitational 
acceleration, r rad = max(g rad ) / g , as a function of effective temperature and gravity. Models 
tend to become numerically unstable when r ra d > 0.9, setting thus limits to our grid. Static 
model atmospheres are increasingly unrealistic so close to the Eddington limit (r rad = 1). 

Finally, we have always adopted a microturbulent velocity, £ t = lOkm/s, which is a 
value typically found in spectral analyses of O stars. These analyses have been performed 
with static model atmospheres as well as with wind codes. Bouret et al. (2003) found that 
Tlusty and CMFGEN models require similar microturbulent velocities to match the UV 
metal lines in SMC O dwarfs, showing therefore that microturbulence is not an artifact 
of static model atmospheres but that it indicates the presence of additional non-thermal 
motions in the atmosphere. We have included the resulting turbulent pressure term in the 
hydrostatic equilibrium equation. 



5.1. Output Products and Availability 

The model atmospheres are available at http://tlusty.gsfc.nas.gov. Each model 
is characterized by a unique filename describing the parameters of the model, for example 
G35000g400vl0. The first letter indicates the composition (see Table 2), followed by the 
effective temperature, the gravity and the turbulent velocity. Each model comes as a set of 
six files, with an identical filename's root but a different extension. A complete description 
of the files' content and format can be found in Tlusty User's guide (see Tlusty web site); 
for reference, we describe them only very briefly here: 

model . 5: General input data; 
model. nst: Optional keywords; 

model . 7: Model atmosphere: Temperature, electron density, total density and NLTE popu- 
lations as function of depth; 

model. 11: Model atmosphere summary; 

model . 12: Model atmosphere: Similar to model . 7, but NLTE populations are replaced by 
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6-factors; 

model. flux: Model flux distribution from the soft X-ray to the far-infrared given as the 
Eddington flux 4 H v [in ergs -1 cm -2 Hz -1 ] as function of frequency. 

The model fluxes are provided at all frequency points included in the calculations (about 
180 000 to 200 000 points with an irregular sampling). Additionally, we have computed 
detailed emergent spectra with Synspec, version 45, in the far-UV (AA900-2000 A) and in the 
optical (AA3000-7500 A). Filename extensions are model. uv. 7, model. vis. 7, respectively, 
and * . 17 for the continuum spectra. Only the optical spectra have been computed for 
the metal-free model atmospheres (model.hhe.7, model .hhe . 17). Additional spectra in 
other wavelength ranges, with altered chemical compositions, or with different values of the 
microturbulent velocity can be readily computed using SYNSPEC. This requires three input 
files, model . 5 , model . 7 , and model . nst, and the necessary atomic data files (model atoms 
and the relevant linelist). 

5.2. Interpolation in the Grid 

A grid of model atmospheres is a good starting point for stellar spectroscopic analyses. 
We can use it to study trends in spectral features with stellar parameters. However, the 
grid sampling is limited by the total number of model atmospheres that can be computed 
in some finite period of time. Therefore, we often have to resort to interpolate in the grid 
for analyzing individual stars. Since radiative transfer and the complete model atmosphere 
problem are highly nonlinear problems, interpolation will result in some errors. We need 
to estimate these errors and find the best way to interpolate in the grid. Here, we consider 
two possibilities: (i) we interpolate the model atmosphere (i.e. the temperature and density 
stratification, as well as the NLTE populations), and then we recompute a spectrum with 
Synspec, or (ii) we interpolate directly the detailed spectra, weighted by (T e g l /T e ff ) 4 . The 
four models and spectra are otherwise given the same weight. 

The Ostar2002 grid is already finely sampled. We decided therefore to estimate the 
interpolation errors by comparing the "exact" UV and optical spectra with spectra interpo- 
lated following the two outlined procedures. We have performed this test on three models, 
S32500g350vl0, L42500g400vl0, and G50000g425vl0. We use the four neighboring models 
(±2 500 K; ±0.25 dex) to obtain the interpolated spectra. In practice, the interpolated model 
will always be closer to a grid point than in these test cases which might therefore be viewed 



4 The flux at the stellar surface is F v — 4tt H v . 
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as worst cases. 

The two interpolations provide spectra almost identical to the "exact" spectra. The full, 
detailed comparison requires too much space to appear here, but the figures are available 
at the Tlusty web site. Minor differences are still found. For example, the spectrum 
interpolation results in a few lines being predicted too strong, especially in the S32500g350vl0 
case. These lines are much stronger in the coolest model spectrum due to lower ionization, 
and this demonstrates the limit of linear spectrum interpolation. On the other hand, we 
found that the model atmosphere interpolation results in too strong N IIIA4640 emissions 
in the L42500g400vl0 models. These emission lines are very sensitive to the exact 6-factor 
ratios, and small differences related to interpolations may result in larger spectral changes. 
Based on these limited tests, we may however conclude that the two interpolation procedures 
appear both reasonably safe because of the sufficiently dense sampling of the Ostar2002 
grid. Nevertheless, we feel that model atmosphere interpolation followed by recomputing 
the spectrum with Synspec remains the safer choice. 

The same helium abundance, solar-scaled compositions, and microturbulent velocity 
have been adopted for all Ostar2002 models. Synspec spectra can nevertheless be re- 
computed with different abundances or microturbulence, using a model atmosphere with the 
appropriate stellar parameters. As long as these new values are not changing drastically 
the overall opacity and thus potentially changing the atmospheric structure, this procedure 
should result in realistic spectra. 

6. Representative Results 

In this section we show a few representative results. Our aim here is not to provide 
a comprehensive display of various models, but instead to show basic trends of important 
model parameters (temperature structure, emergent spectra) with effective temperature, 
surface gravity, and metallicity. 

Fig. 4 shows the local temperature for a constant surface gravity (log g = 4), for three 
values of T e g (30 000, 40 000 and 50 000 K), and for all values of metallicity. The temperature 
distribution nicely illustrates the basic features of line blanketing, namely the backwarming - 
line blanketing leads to a heating of an atmosphere roughly between Rosseland optical depth 
0.01 and 1; and the surface cooling. The zero- and very low- metallicity models exhibit a 
temperature rise at the surface, a typical NLTE effect (an indirect heating effect of hydrogen 
Lyman and Balmer lines), discovered and explained already in early days of NLTE model 
atmospheres by Auer & Mihalas (1969a). This effect competes with surface cooling due to 
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metal lines. At 50 000 K, both effects nearly cancel at metallicity 1/50 of solar, while at 
lower T eS the rise is wiped out at higher metallicities (around 1/10 of solar). Interestingly, 
the temperature curves for all metallicities cross in a narrow range of optical depths; this 
crossing depth decreases with decreasing T eS (being about 0.025 at 50 000 K; 0.02 at 40 000 K, 
and 0.002 at 30 000 K). 

Figures 5-9 display the ionization fractions of all explicit species, along with LTE frac- 
tions (for a direct comparison, the LTE fractions are computed by the Saha formula for 
temperature and density distributions of the NLTE models - that is, we did not recompute 
self-consistent LTE model structures). The behavior of individual species is different, but 
we can see a general trend: the fractions of dominant ions are usually not sensitive to NLTE 
effects (that is, they are close to unity in both LTE and NLTE). The fractions of ions with 
lower ionization degrees are typically lower in NLTE, because these ions are overionized by 
the strong radiation field that originates in deep, very hot layers (notice that the recombina- 
tion rate is given through the local temperature). Conversely, higher ionization degrees are 
typically larger in NLTE. This ionization shift is one of the most important NLTE effects in 
hot star atmospheres. 

Next, we study a sensitivity of emergent flux on T c s, log g, and metallicity. For simplicity, 
we display spectra computed directly by Tlusty (files model. flux described earlier), not 
detailed spectra computed by Synspec. For clarity, the spectra are smoothed over 500 
frequency points, which roughly simulates a 10 A resolution of Kurucz (1993) model spectra. 

Fig. 10 shows a sequence of model UV spectra with logg = 4 and solar metallicity, 
for various effective temperatures. The Lyman jump gradually weakens with increasing 
effective temperature, and essentially disappears at 50 000 K. Similarly, one can easily observe 
gradual weakening and ultimate disappearance of the Si IVA1400, C IIIA1176, 977, Lya, 
etc., features. Fig. 11 shows the effect of surface gravity for a given effective temperature 
(40 000 K) and metallicity (solar). The most striking feature is the magnitude of the Lyman 
jump, which is very strong at logg = 4.5, and almost disappears at log g = 3.5. Similarly, 
the line strengths vary significantly. Next, we show the effect of metallicity. Fig. 12 displays 
models for T cS = 40000 K, \ogg = 4, and metallicities Z/Z Q = 2, 1, 0.5, 0.3, and 0.1. 
In the low-metallicity models, the Lyman continuum flux and the near-UV continuum flux 
(A > 1700 A) is lower than that for higher metallicities. This relation is reversed in the 
region crowded with metal lines, 900 < A < 1700 A. 

Finally, we compare our model fluxes to Kurucz (1993) models. Because of different 
averaging procedures used in both sets of model spectra, the comparison is not exact, nev- 
ertheless shows the general effect. We compare three models for (T e g, log g) pairs equal to 
(40 000, 4.5), (35 000, 4.0), and (30 000, 4.0), and for solar metallicity - Fig.13. The agree- 
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ment between our and Kurucz models is very good in the optical spectrum for all three 
models (lower panel). For T e g = 30 000 K, the agreement between Ostar2002 and Kurucz 
fluxes is rather good in the UV spectrum (A > 900 A), while for the hotter models our NLTE 
flux is lower. In contrast, NLTE models predict higher flux than Kurucz model in the Lyman 
continuum, although the effect is relatively small at 40 000 K (see also §7). 

This exercise shows that using Kurucz model is still a reasonable choice for analyzing 
low-resolution UV and optical spectra of O stars. It should be realized, however, that the 
full strength of NLTE models lie in their ability to provide reliable high-resolution spectra 
and profiles of individual lines that can be used for determining the stellar parameters and 
the surface chemical composition. We will study various systematic effects and differences 
between LTE and NLTE models, and their influence on deduced stellar parameters, in a 
future paper. 

7. Bolometric Corrections and Ionizing Fluxes 
7.1. Bolometric Corrections 

Bolometric corrections (Table 3) are calculated using the following expression: 

BC = m hol - V = (-2.51ogF bo i - 11.487) - (-2.5 \ogF v - 21.100) (1) 

where F bo i is the bolometric flux and F v is the flux through the Johnson V filter, com- 
puted from the model atmospheres' flux distribution. The bolometric flux is computed by 
trapezoidal integration over the complete frequency range, while we have used an IDL ver- 
sion of the program UBVBUSER (Kurucz 1993) kindly made available by Wayne Landsman 
to obtain the V magnitudes. The first constant is defined by assuming a solar constant, 
f® ol = 1.371 10 6 erg s _1 cm~ 2 , a solar visual magnitude, V & = —26.76, and a solar bolomet- 
ric correction, BC Q = —0.07, while the second constant defines the zero point of the V 
magnitude scale (Bessel et al. 1998). 

Fig. 14 illustrates the major dependence of the bolometric correction with the effective 
temperature. Over the parameter space of the grid, a multi-linear regression yields the 
relation: 

BC = 27 A3 - 6.78 x log(T cff ) + 0.06 x Z/Z Q (2) 

An inspection of Table 3 shows that the dependence with gravity is generally even smaller 
than the metallicity term. We have therefore neglected the gravity dependence and Eq. 2 has 
been established with values from log g = 4.0 models. The typical accuracy of bolometric 
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corrections derived from Eq. 2 is ±0.05 mag. Larger differences up to 0.1 mag between the 
models and Eq. 2 values are found for the coolest models. 

We compare our bolometric corrections to values calculated for Kurucz (1993) models 
and extracted from Bessel et al. (1998). We use the same zero point for bolometric correc- 
tions. We find almost identical values for solar composition models with logg = 4.0, and 
4.5. The largest differences reach 0.03 mag at T cS = 47500 K, and 0.05 mag at T cS = 27500 K, 
adding thus to the good agreement found in the previous section between spectral energy 
distributions (Fig. 13). 

Our relation (Eq. 2) is close to the Vacca, et al. (1996) calibration, which is based on a 
number of diverse NLTE analyses. Their bolometric corrections are systematically smaller 
by 0.12 mag. Our relation therefore yield stellar luminosities 12% fainter than previously 
estimated from Vacca et al. (1996). 

7.2. Ionizing Fluxes 

We have calculated the ionizing fluxes in the hydrogen Lyman and in the He IA504 
continua. They are listed in Tables 4 and 5. The ionizing fluxes, go an d gi, are given 
as logarithms of the number of photons in these two continua, per second and per square 
centimeter at the stellar surface. We do not provide the ionizing flux q 2 in the He II Lyman 
continuum. Gabler et al. (1991) indeed showed that the He II continuum is formed in a stellar 
wind (due to high opacity in this continuum), and that spherical extension and shocks are 
necessary model ingredients to predict fluxes that are in rough agreement with observed 
fluxes. However, the contribution of the He II Lyman continuum to g and gi remains very 
small. Any error in g 2 results therefore in unsignificant changes of go and gi. The latter 
fluxes are mostly sensitive to a correct treatment of NLTE metal line blanketing, which in 
turn is crucial to predict the hydrogen and helium ionization structure and the photospheric 
temperature stratification reliably. 

The dependence of the ionization fluxes with T c fj, logg, and Z/Z Q , is displayed in 
Fig. 15. Lower gravity models have larger ionizing fluxes due to higher hydrogen and helium 
ionization (thus the Lyman and A504 jumps are smaller). Model atmospheres with low 
metallicity have lower Lyman continuum fluxes. The behavior of the He I continuum fluxes 
is more complicated: at low T e fj, metal-poor model atmospheres have markedly lower qi 
fluxes, while at T c s higher than about 33 000 K, they have slightly larger gi fluxes. Lower 
metal abundances reduce the blanketing effect, and thus result in lower temperatures in the 
continuum-forming layers and lower fluxes. The higher qi fluxes in hot, metal-poor models 
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come from smaller metal line opacities in this spectral range. 

We have compared our ionizing fluxes to other recent works. Vacca et al. (1996) used 
Kurucz (1993) LTE model fluxes to compute q and q\. They list 43 values for parameters 
covering the domain of O stars. We interpolate bilinearly in logT cff and logg in Tables 4 
and 5 and compare to Vacca et al.'s values. The qo values agree very well, with no systematic 
differences; we still found differences up to 0.1-0.2 dex for some parameters (mostly the coolest 
models). On the other hand, Vacca et al.'s q± values are on average 0.25 dex smaller than 
ours; individual differences range typically from 0.1 to 0.5 dex, increasing towards cooler 
effective temperatures. These differences are quite significant for nebular studies. 

We have also compared our results to values derived from NLTE, line-blanketed, ex- 
panding model atmospheres, calculated with program wm-basic (Pauldrach et al. 2001). 
The go values agree also very well, although our values are slightly larger. The values agree 
within 0.1 dex, with the exception of the two models, D-30 and S-30, at T e s= 30 000 K, where 
wm-basic models predict ionizing fluxes smaller by 0.57 and 0.24 dex, respectively. 

8. Model Sensitivity to Atomic Model Atoms 

We have aimed at incorporating all important opacity sources and all relevant atomic 
processes that determine the excitation and ionization balance in the model atmospheres. 
To investigate the sensitivity of the model atmospheres to our handling of atomic data, we 
present here a few test cases where we adopted different model atoms. We consider here 
two potential issues: (i) the size of the iron model atoms and the representation of the line 
opacity; (ii) the omission of highly-excited levels in light elements (mostly, CNO). 

Kurucz (1994) lists thousands of energy levels for every iron ion. How should we group 
them? Our model atoms group them into 30 to 50 superlevels per ion. On the other hand, 
Deetjen et al. (1999) have grouped all individual levels into 7 to 9 NLTE superlevels for 
each iron ion. We have compared the effect of different Fe model data on the resulting 
model atmospheres, with a test case, T e s = 35 000 K, log g = 4.0, solar composition, and 
£t = 10 km/s. The reference model is extracted from the Ostar2002 grid. The first three 
models deal with the representation of line opacity: OS with larger frequency step, ODF, 
and line strength selection (Fe and Ni lines are dynamically selected depending on a strength 
criterion, see §4.10). In the case of the ODF model (Gb), all the Fe and Ni lines are used to 
set up the ODFs. In the second set of models, we have built different Fe IV and Fe v model 
atoms (they are the dominant ions, see Fig. 9). All individual levels are included in models 
Gd and Ge, but grouped in more (74&69 vs. 43&42) or less (33&19) superlevels. In the last 
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model, all the Fe and Ni levels above the ionization limits are omitted. Table 6 summarizes 
the properties of the different models. 

Fig. 16 illustrates the changes in the model temperature atmosphere structure due 
to different ways of including the iron data. This is a condensed way to assess differences 
between model atmospheres. We may indeed expect that model spectra change with changes 
in the temperature structure (this does not exclude of course that some line profiles might 
be different even with very close temperature structures, see below). First, differences are 
larger with ODF's than when adopting a rather large sampling step. We interpret this 
as reflecting incorrect blends between ODFs and important lines of light elements. The 
larger sampling step is probably appropriate in many cases (AT 100 K), and reduces the 
computing time by a factor of about three (the reduction in the number of frequency points). 
The largest difference arising due to the line strength selection criterion is seen at depth, 
because of a large change in the number of selected Fe VI lines. However, this does not 
affect the predicted spectrum. Differences in the layers where the spectrum is formed are 
small (AT < 100 K) showing that most of the line-blanketing effect is due to the strongest 
10 to 10 5 lines rather than to millions of weaker lines. The lower panels of Fig. 16 show the 
results related to the model atoms. First, we see that increasing the number of superlevels 
results in small changes, while the model with less superlevels shows a larger systematic 
temperature difference. Although limited in scope (one set of stellar parameters), this test 
thus supports a claim that we have adopted Fe model atoms of reasonable sizes for our grid 
of model atmospheres. Finally, the levels above the ionization limits have a minimal effect 
on the model atmosphere structure. 

In a second test, we have investigated the omission of highly-excited levels of light 
elements. Sect. 4 details the adopted model atoms. In some cases, levels with main quantum 
number, n > 7, have been neglected (e.g., C III and N IV). We have built a new C III 
model atom which includes all levels below the ionization limit that are considered in OP 
calculations. The new model atom includes the first 34 levels as individual explicit levels 
and group higher levels into 12 superlevels. We have recomputed one model atmosphere, 
S35000g325vl0, assuming the extended C III model atom. The atmospheric structure is 
changed very little (AT < 8K), and this omission is therefore mostly of no consequence. 

However, the total recombination rate, from C +3 to C +2 is significantly increased, by a 
factor 2 to 3. Recombination lines in the optical, e.g. C IIIA4650, show thus significantly 
stronger emission in the new model. We note here in passing that these "Of-type" lines 
are predicted in emission by plane-parallel, hydrostatic, NLTE model atmospheres at low 
gravities (Mihalas & Hummer 1973). These are NLTE recombination lines and emission is 
neither a consequence of a stellar wind nor of an extended atmosphere. We caution therefore 
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against using the current grid for quantitative analysis of weak absorption or emission fea- 
tures in the optical spectrum without a careful investigation case by case. A future update 
is expected. 

9. Discussion and Conclusions 

We have constructed a comprehensive grid of metal line-blanketed, NLTE, plane-parallel, 
hydrostatic model atmospheres for the basic parameters appropriate to O-type stars. The 
Ostar2002 grid considers 12 values of effective temperatures 27 500 K <T cS < 55 000 K, 
with 2 500K steps, 8 surface gravities, 3.0 < logg < 4.75 with 0.25 dex steps, and 10 chemi- 
cal compositions, from metal-rich relative to the Sun to metal-free. The lower limit of logg 
for a given T e g is actually set by an approximate location of the Eddington limit. The se- 
lected chemical compositions have been chosen to cover a number of typical environments of 
massive stars: the galactic center, the Magellanic Clouds, Blue Compact Dwarf galaxies like 
IZwl8, and galaxies at high redshifts. The paper contains a description of the Ostar2002 
grid and some illustrative examples and comparisons. The complete Ostar2002 grid is 
available at our website at http://tlusty.gsfc.nasa.gov. 

We are aware that the presented grid cannot represent the last word in O star model 
stellar atmospheres. However, we intended to provide a more or less definitive grid of model 
atmospheres in the context of one-dimensional, plane-parallel, homogeneous, hydrostatic 
models in radiative equilibrium. We have attempted to take into account essentially all im- 
portant opacity sources (lines and continua) of all astrophysically important ions. Likewise, 
we have attempted to consider all relevant atomic processes that determine the excitation 
and ionization balance of all such atoms and ions. In other words, since the computational 
effort behind constructing the above described 690 models is enormous, we intended to com- 
pute as complete and realistic models as possible, in order to avoid a need to recompute a 
grid in a near future should some neglected mechanism or an opacity source is found to be 
important. 

Although we spent all effort to make sure that the treatment of atomic physics and opac- 
ities is as complete and accurate as possible, there are still several points we are aware of that 
were crudely approximated. Those approximations were necessitated not by shortcomings 
of our modeling scheme or a lack of adequate computer resources, but by the present lack 
of sufficient atomic data. These approximations include (roughly in the order of estimated 
importance): 

• A crude and approximate treatment of collision rates for forbidden transitions of Fe 
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and Ni. As discussed in § 4, since no approximate formula or a simple prescription is available, 
we have somewhat arbitrarily adopted the Eissner-Seaton formula with 7(T) = 0.05. The 
IRON project provides, and will keep providing in future, theoretical calculations for selected 
transitions for various ions of iron; we plan to implement all the newly available collision 
rates in future calculations in order to check their influence explicitly. 

• A neglect of very high energy levels of light metals. The Opacity Project, which is the 
basic source of our atomic data, typically considers energy levels only up to a certain principal 
quantum number (usually 10). The higher levels are neglected. This does not usually cause 
any significant errors in the excitation balance, but in some cases recombinations to these 
neglected levels may represent a non-negligible contribution to the total recombination rate 
- the global ionization balance may thus be somewhat in error. 

• Inaccuracies in, or a lack of, computed photoionization cross-sections. Although the 
Opacity Project provided a large number of cross-sections, many are still missing. We have 
approximated all cross-sections for which we did not find any published data by hydro- 
genic cross-sections. However, since we have adopted accurate data for the most important 
transitions, we do not believe that this may be a source of significant errors. 

• We did not consider charge exchange reactions. In this case, some necessary atomic 
data are available, but in view of high temperature of the atmospheres of O-type stars, 
we do not expect the charge exchange reactions to have any significant contribution to the 
computed NLTE excitation/ionization balance. 

For completeness, we repeat the list the most important physical approximations of our 
model grid: 

• Plane-parallel atmosphere in hydrostatic equilibrium. This is perhaps the most im- 
portant and restrictive assumption. Because of its importance, we have discussed it at length 
in §2.1. Essentially, we argued that while the wind features cannot obviously be described 
by means of our models, the rest of the UV and optical spectra is well described, unless the 
wind is truly extreme. In any case, the Ostar2002 grid should not be used for analyzing 
Luminous Blue Variables or Wolf-Rayet stars, where already the continuum is formed in the 
wind; and the grid should be used with caution in the case of Of supergiants. 

• Homogeneity (i.e., a 1-D atmosphere). This means that any surface inhomogeneities 
are neglected. However, our understanding of inhomogeneities in hot star atmospheres is 
very poor, and at present there is little what can be done on a quantitative level. 

• Radiative equilibrium. Our approach neglects any source of mechanical energy dissi- 
pation, shocks, non-thermal X-rays, and related phenomena. Again, there is no approach at 
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present that would address those issues from the first principles, and therefore any modeling 
of such phenomena is necessarily extremely crude. 

Despite of the long list of approximations and uncertainties, we believe that the Os- 
tar.2002 grid represents a definite improvement over previous grids of O star model atmo- 
spheres. We believe that we have constructed a more or less definitive grid of models in 
the context of one-dimensional, plane-parallel, homogeneous, geometry in a hydrostatic and 
radiative equilibrium, with all important opacity sources (metal lines) included, and without 
any unnecessary numerical approximations. We hope that the OSTAR2002 grid of models 
will be useful for a number of years to come, until a new-generation models with multi- 
dimensional, self-consistent radiation magneto-hydrodynamic description will be developed. 

We are indebted to Sally Heap for her continuing encouragement and motivation to 
compute this grid, and for suggesting a number of various improvements of the modeling 
strategy. We acknowledge helpful discussions and collaboration with many colleagues, most 
importantly with John Hillier, David Hummer, Dimitri Mihalas, Klaus Werner, Stefan Drei- 
zler, Stefan Haas, and the late Keith Smith. On the computational side, we would like to 
thank most sincerely Keith Feggans for his invaluable assistance, and for his patience to 
answer our many calls to help whenever any computer-related problem occurred, day or 
night, during our work on computing model atmospheres at NASA/GSFC. Last but not 
least, we are greatly indebted to all the physicists who provided, and are still providing, 
all the necessary atomic data, in particular the OP/IP team, Bob Kurucz, and the NIST 
team. The present sophistication of current model atmospheres would simply not been 
achievable without all their work during many years. This work was supported by NASA 
grant NAG5-10895 (FUSE B134 program), NASA NRA-01-01-ADP-099 grant, and by sev- 
eral grants (GO 7437, AR 7985, GO 9054) from the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS5-26555. 
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Fig. 1. — Selected Ostar2002 grid points in the \ogg vs. T cS plane. Geneva evolutionary 
tracks (Schaller et al. 1992; Charbonnel et al. 1993) are shown for solar and 1/5 solar metal- 
licities (full and dashed lines, respectively). The tracks correspond to models with initial 
masses of 120, 85, 60, 40, 25, 20, 15, 12 M Q from left to right, respectively. 
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Fig. 2. — Radiative acceleration (thick line) for a solar composition model atmosphere, T c g 
= 40 000, log g = 4.0. At low optical depths, (r ross < 10~ 6 ), the total radiative acceleration 
(dashed line) increases and becomes larger than the gravitational acceleration (thin line). We 
apply a numerical cut-off to the radiative acceleration to ensure convergence. The assumption 
of hydrostatic equilibrium breaks down in these superficial layers. The arrow indicates our 
definition of stability against radiation pressure (r rac j < 1). 
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Fig. 3. — Isolines of radiative acceleration relative to gravitational acceleration, r ra d = gr&d/g, 
for solar composition model atmospheres, as a function of T cff and log p. r rad values are listed 
right to the isolines. The dashed line shows an estimate of the Eddington limit obtained by 
extrapolation. 
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Fig. 4.— Model atmospheres for T eff = 30 000, 40 000 and 50 000 K (top to bottom panels), 
logg = 4.0, and for various metallicities. At low optical depths (r ross < 10~ 3 ), the top 
curves is for a H-He model, and the temperature is progressively lower when increasing the 
metallicity, while reverse is true at deeper layers (r ross > 10~ 2 ). 
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Fig. 5. — Ionization fractions of helium and carbon in three model atmospheres, T e g= 
30 000, 40 000 and 50 000 K (top to bottom panels), \ogg = 4.0, and solar composition. LTE 
ionization is shown with dashed lines. 
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Fig. 6. — Ionization fractions of nitrogen and oxygen in three model atmospheres, T e g= 
30 000, 40 000 and 50 000 K (top to bottom panels), \ogg = 4.0, and solar composition. LTE 
ionization is shown with dashed lines. 
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Fig. 7. — Ionization fractions of neon and silicon in three model atmospheres, T eS = 30 000, 
40 000 and 50 000 K (top to bottom panels), log gi = 4.0, and solar composition. LTE ioniza- 
tion is shown with dashed lines. 
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Fig. 8. — Ionization fractions of phosphorus and sulfur in three model atmospheres, T e g= 
30 000, 40 000 and 50 000 K (top to bottom panels), \ogg = 4.0, and solar composition. LTE 
ionization is shown with dashed lines. 
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Fig. 9. — Ionization fractions of iron and nickel in three model atmospheres, T eS = 30 000, 
40 000 and 50 000 K (top to bottom panels), \ogg = 4.0, and solar composition. LTE ioniza- 
tion is shown with dashed lines. 
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Fig. 10. — Predicted flux for model atmospheres for log g = 4.0, solar composition, and 
various effective temperatures - from top to bottom T eS = 55, 50, 45, 40, 35, and 30 kK. 
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Fig. 11. — Predicted flux for model atmospheres for T cS = 40 000 K, solar composition, and 
various surface gravities - from top to bottom at the Balmer continuum (A > 912 A) log g = 
4.5, 4.25, 4.0, 3.75, and 3.5; the lowest-gravity model is drawn by bold line. Notice that the 
dependence of flux on logg is reversed in the Lyman continuum (A < 912 A). 
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Fig. 12. — Predicted flux for model atmospheres for T cS = 40 kK, logg = 4.0, and various 
metallicities - Z/Z Q = 2, 1, 1/2, 1/5, 1/10. The flux for the lowest-metallicity model is drawn 
by a bold line. The flux generally decreases with decreasing metallicity for A < 900 A, and 
A > 1700 A, while the relation is reversed for 900 < A < 1700 A. 
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Fig. 13. — Predicted flux for model atmospheres for three model atmospheres with (T eff , log g) 
equal to (40000 K, 4.5); (35000 K, 4.0), and (30000 K, 4.0) - bold lines; compared to Kurucz 
models with the same parameters - thin histograms. Different panels show different spectrum 
regions. 
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Fig. 14. — Bolometric corrections vs. effective temperature. The top panel displays the 
relation for \ogg = 4.0, solar composition (full), 1/10 (dashed) and 1/1000 solar (dotted). 
The middle panel expands on the effect of metallicity on BC, displaying differences relative 
to solar models (log g = 4.0). The bottom panel shows the effect of different gravities (full: 
solar; dashed: 1/10 solar). 
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Fig. 15. — Ionization fluxes [s _1 cm~ 2 ] vs. effective temperature. Top: log of number of 
Lyman continuum photons; bottom: log of number of He IA 504 photons. In each panel, there 
are 4 families of curves from bottom to top with decreasing gravities, log g = 4.5, 4.0, 3.5, 3.0. 
Results for different chemical compositions are shown with different types of lines: solar (full), 
1/10 solar (dashed), and metal-free (dotted). 
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Fig. 16. — Changes in the model temperature structure with respect to the reference model 
depending on the treatment of Fe model atoms and line opacity (see Table 6 for model keys). 
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Table 1. Atomic data included in the NLTE model atmospheres. 



Ion (Super)Levels Indiv. Levels Lines References 



HI 


9 


80 


172 




HII 


1 


1 






He I 


24 


72 


784 


Martin et al. (1999) 


He II 


20 


20 


190 




He III 


1 


1 






CII 


22 


44 


238 


Fernley et al. (1999) 


cm 


23 


55 


258 


Tully et al. (1990) 


CIV 


25 


55 


330 


Peach et al. (1988) 


C V 


1 


1 






Nil 


26 


71 


373 


Luo k Pradhan (1989) 


NIII 


32 


68 


549 


Fernley et al. (1999) 


NIV 


23 


58 


286 


Tully et al. (1990) 


N V 


16 


55 


330 


Peach et al. (1988) 


N VI 


1 


1 






on 


29 


62 


346 


Burke k Lennon (unpub.) 


III 


29 


116 


846 


Luo k Pradhan (1989) 


OIV 


39 


94 


922 


Fernley et al. (1999) 


V 


40 


89 


610 


Tully et al. (1990) 


VI 


20 


55 


330 


Peach et al. (1988) 


VII 


1 


1 






Ne II 


15 


29 


83 


Butler k Zeippen (unpub. 


Ne III 


14 


20 


32 


Butler k Zeippen (unpub. 


Ne IV 


12 


18 


38 


Burke k Lennon (unpub.) 


Ne V 


1 


1 






Si III 


30 


105 


747 


Butler et al. (1993) 


Si IV 


23 


53 


306 


Taylor (unpub.) 


Si V 


1 


1 







-49 - 



Table 1 — Continued 



Ion (Super)Levels Indiv. Levels Lines References 



P IV 14 14 19 Martin et al. (1999) 

PV 17 43 223 

P VI 1 1 

S III 20 27 57 Nahar & Pradhan (1993) 

S IV 15 15 26 Mendoza et al. (1995) 

S V 12 12 11 Butler et al. (1993) 

S VI 16 34 147 Taylor (unpub.) 

S VII 1 1 

FeUI 50 12 660 1604 934 Kurucz (1994); Nahar (1996) 

FelV 43 13 705 1776 984 Kurucz (1994); Bautista & Pradhan (1997) 

FeV 42 11986 1008 385 Kurucz (1994); Bautista (1996) 

FeVI 32 4 740 475 750 Kurucz (1994) 

Fe VII 1 1 

NiHI 36 11335 1309 729 Kurucz (1994) 

NilV 38 13172 1918 070 Kurucz (1994) 

NiV 48 13184 1971819 Kurucz (1994) 

NiVI 42 13 705 2 259 798 Kurucz (1994) 

Ni VII 1 1 
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Table 2. Key to models chemical compositions. 



Key 


Composition 


Key 


Composition 


C 


2x0 


V 


1/30 x 


G 


1 x 


w 


1/50 x 


L 


1/2 x 


X 


1/100 x 


S 


1/5 x 


Y 


1/1000 x 


T 


1/10 x 


Z 
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Table 3. Bolometric corrections as function of effective temperature, gravity and 
metallicity (10 metallicities from 2 times solar to metal-free models). 



T cS [K] 


log g 


2. 


1. 


0.5 


0.2 


BC 
0.1 


[mag] 
0.033 


0.02 


0.01 


0.001 


0. 


z / ouu 


i nn 
o.uu 


9 ^Q 

— z.oy 


9 P.0 

— z.uz 


9 px, 
— Z.DO 


9 PR 
— Z.Do 


9 71 
— Z. 1 1 


9 71 
— Z. / O 


9 7/1 
— Z. / 4 


9 7X. 
— Z. / O 


9 7P 
— Z. f O 


9 7P. 
— Z. 1 


97^nn 
Z 1 ouu 


1 9^ 
O.ZO 


9 x.c 
— Z.OD 


9 fil 
— Z.D1 


9 px. 
— Z.DO 


9 fiQ 

— z.oy 


9 71 
— Z. 1 1 


9 71 
— Z. / O 


9 7/1 
— Z. / 4 


9 7/1 
— Z. / 4 


9 7X. 
— Z. 1 O 


9 7X. 
— Z. ( 


Z 1 ouu 


i ^n 
o.ou 


9 ^7 
— Z.O / 


9 fi9 

— z.uz 


9 PP 
— Z.DO 


9 7n 

— Z. ( u 


9 79 
— Z. / Z 


9 7/1 
— Z. / 4 


9 7^ 
— Z. / O 


9 7fi 
— Z. / O 


9 7P 
— Z. f O 


9 7P. 
— Z. ( O 


z / ouu 


1 7X. 
O. f O 


9 ^Q 

— z.oy 


9 P.A 
— Z.04 


9 fi8 
— Z.Do 
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Table 3 — Continued 



T cS [K] 


log g 


Z. 


1. 


0.5 


0.2 


BC 

n 1 
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[mag] 
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n no 
0.02 


0.01 


0.001 


U. 
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-3.69 
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-3.89 
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-4.26 
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-4.18 


-4.20 


-4.22 


-4.24 


-4.25 
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-4.64 


-4.66 


-4.68 
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-4.71 
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Table 4. Ionizing fluxes in the H I Lyman continuum as function of effective temperature, 
gravity and metallicity (10 metallicities from 2 times solar to metal-free models). 
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Table 5. Ionizing fluxes in the He I continuum as function of effective temperature, 
gravity and metallicity (10 metallicities from 2 times solar to metal-free models). 
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Table 6. Assumptions and characteristics for the test case models. 



Model 


Characteristics 


Levels 


Freqs 


Fe Lines 


G35000g400vl0 


Reference model (OS, 0.75) 


907 


184136 


1 176 853 


Ga35 


OS, 30 Doppler widths 


907 


65 398 


1 176 821 


Gb35 


ODF 


907 


40 759 


7459 416 


Gc35 


Fewer selected Fe/Ni lines 


907 


184136 


88 009 


Gd35 


More Fe iv-v superlevels 


965 


65 493 


1221466 


Ge35 


Less Fe iv-v superlevels 


874 


65 366 


952 051 


Gf35 


No levels above ioniz. limit 


804 


184226 


1 168 764 



